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Stochastic Background of Relic Scalar Gravitational Waves tuned by 
Extended Gravity 

Mariafelicia De Laurentis a , Salvatore Capozziello a 

a Dipartimento di Scienze Fisiche, Universita di Napoli " Federico II", 

INFN Sez. di Napoli, Compl. Univ. di Monte S. Angelo, Edificio G, Via Cinthia, 1-80126, Napoli, Italy 

A stochastic background of relic gravitational waves is achieved by the so called adiabatically-amplified zero- 
point fluctuations process derived from early inflation. It provides a distinctive spectrum of relic gravitational 
waves. In the framework of scalar-tensor gravity, we discuss the scalar modes of gravitational waves and the 
primordial production of this scalar component which is generated beside tensorial one. Then we analyze different 
viable /(i?)-gravities towards the Solar System tests and stochastic gravitational waves background. The aim is 
to achive experimental bounds for the theory at local and cosmological scales in order to select models capable 
of addressing the accelerating cosmological expansion without cosmological constant but evading the weak field 
constraints. It is demonstrated that viable /(-R)-gravities under consideration not only satisfy the local tests, but 
additionally, pass the PPN-and stochastic gravitational waves bounds for large classes of parameters. 



1. Introduction 

The idea that General Relativity (GR) should 
be extended or corrected at large scales (infrared 
limit) or at high energies (ultraviolet limit) is sug- 
gested by several theoretical and observational is- 
sues [ij. Quantum field theory in curved space- 
times, as well as the low-energy limit of String/M 
theory, both imply semi-classical effective actions 
containing higher-order curvature invariants or 
scalar-tensor terms. In addition, GR has been 
definitely tested only at Solar System scales while 
it may show several shortcomings if checked at 
higher energies or larger scales. Besides, the Solar 
System experiments are, up to now, not conclu- 
sive to state that the only viable theory of gravity 
is GR: for example, the limits on PPN parameters 
should be greatly improved to fully remove degen- 
eracies. Of course, modifying the gravitational 
action asks for several fundamental challenges. 
These models can exhibit instabilities or ghost- 
like behavior, while, on the other hand, they have 
to be matched with observations and experiments 
in the appropriate low energy limit, there exist 
cosmological solutions that give the accelerated 
expansion of the universe at late times. In ad- 
dition, it has been discovered that some stability 



conditions can lead to avoid ghost and tachyon 
solutions. Furthermore there exist viable f(R) 
models which satisfy both background cosmolog- 
ical constraints and stability conditions and re- 
sults have been achieved in order to place con- 
straints on f(R) cosmological models by CMBR 
anisotropies and galaxy power spectrum. More- 
over, some of such viable models lead to the uni- 
fication of early-time inflation with late-time ac- 
celeration. On the other hand, by considering 
/(i?)-gravity in the low energy limit, it is pos- 
sible to obtain corrected gravitational potentials 
capable of explaining the flat rotation curves of 
spiral galaxies or the dynamics of galaxy clus- 
ters without considering huge amounts of dark 
matter. Furthermore, several authors have dealt 
with the weak field limit of fourth order gravity, 
in particular considering the PPN limit and the 
spherically symmetric solutions. This great deal 
of work needs an essential issue to be pursued: 
we need to compare experiments and probes at 
local scales (e.g. Solar System) with experiments 
and probes at large scales (Galaxy, extragalac- 
tic scales, cosmology) in order to achieve self- 
consistent f(R) models. In order to constrain 
further viable /(i?)-models, one could take into 
account also the stochastic background of gravi- 
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tational waves (GW) which, together with cosmic 
microwave background radiation (CMBR), would 
carry a huge amount of information on the early 
stages of the Universe evolution. In fact, if de- 
tected, such a background could constitute a fur- 
ther probe for these theories at very high red-shift 
[9]. On the other hand, a key role for the produc- 
tion and the detection of the relic gravitational 
radiation background is played by the adopted 
theory of gravity . This means that the effective 
theory of gravity should be probed at zero, inter- 
mediate and high redshifts to be consistent at all 
scales and not simply extrapolated up to the last 
scattering surface, as in the case of GR. 

The aim of this report is to discuss the PPN 
Solar-System constraints and the GW stochastic 
background considering some recently proposed 
f(R) gravity models |2I3| which satisfy both 
cosmological and stability conditions mentioned 
above. Using the definition of PPN-parameters 
7 and in terms of /(i?)-models and the defini- 
tion of scalar GWs [4], we compare and discuss 
if it is possible to search for parameter ranges of 
/(i?)-models working at Solar System and GW 
stochastic background scale [5]. 



2. Extended Gravity 

A action for scalar-tensor gravity in a Brans- 
Dicke-like form which can be adopted also for 
/(i?)-gravity once a suitable scalar field is defined, 

is m 



S = Id x^/—g 



_ v - W(cp) + Cr, 

(1) 



By varying the action (!]) with respect to g M „, we 
obtain the field equations 



Gut, 



i-wGrpjm) u(ip) ( 1 a/3 \ 



(2) 



while the variation with respect to p> gives the 
Klein - Gordon equation 



Dp = 



1 



2u(ip) + 3 

+2W(p) + pW'(p) + ^M.g^ W . v ]. 

(3) 

We are assuming physical units G = 1 , c = 1 and 
h = 1. is the matter stress-energy tensor 

and G is a dimensional, strictly positive, gravita- 
tional coupling constant |4I7| . The Newton con- 
stant is replaced by the effective coupling 



G, 



eff 



1 



(4) 



which is, in general, different from G. GR is re- 
covered for 



(5) 



3. Gravitational waves from Extended 
Gravity 

In order to study gravitational waves, we as- 
sume first-order, small perturbations in vacuum 
(T / i™ ) = 0). This means 



and 



g^v = Viiv +h^ v , tp = tpo + Sp (6) 



W ~ -aStp 2 =>W ~ aSp (7) 

for the self-interacting, scalar-field potential. 
These assumptions allow to derive the "lin- 
earized" curvature invariants R^ upa , R^v and R 
which correspond to R^upa , R^u and R, and then 
the linearized field equations |4I8| 



(8) 



R P u - §Vn» = -d M d v ® + 
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In particular, the transverse-traceless (TT) gauge 
(see |8j) can be generalized to scalar-tensor grav- 
ity obtaining the total perturbation of a GW in- 
coming in the z+ direction in this gauge as 



+$(t-*) e W 



(10) 



The term A+(t - z)e$ + A x (t - z)e { t _ l l' describes 
the two standard (i.e. tensorial) polarizations of 
a gravitational wave arising from GR in the TT 
gauge [8], while the term $(t — z)e|fj is the ex- 
tension of the TT gauge to the scalar case. This 
means that, in scalar-tensor gravity, the scalar 
field generates a third component for the tensor 
polarization of GWs. This is because three differ- 
ent degrees of freedom are present of |4|10j , while 
only two are present in standard General Relativ- 
ity. 

4. Stochastic background of relic scalar 
GWs 

Then, for a purely scalar gravitational wave, 
the metric perturbation is |4|10|9j 



(x) 



(11) 



The stochastic background of scalar gravitational 
waves can be described in terms of the scalar field 
$ and characterized by a dimensionless spectrum 
(see the analogous definition for tensor modes in 
01) 



i dp 



sgw 



Pc din/ 



where 



mi 

8irG 



(12) 



(13) 



is the (present) critical energy density of the 
universe, Hq is the Hubble parameter today, 
and dp sgw is the energy density of the scalar 
gravitational radiation in the frequency interval 
(/, / + d/). We are now using standard units. 
Now it is possible to write an expression for the 
energy density of the stochastic scalar relic gravi- 
tons background in the angular frequency interval 



(u>, lo + duj) as 



dpsgw = 2fiw , , , 

hH 2 s H 2 doj _ hH 2 s H 2 df 



47T 2 C 3 U> 



47T 2 C 3 / 



(14) 



where /, as above, is the frequency in standard 
comoving time. Eq. (fl4|l can be rewritten in 
terms of the critical and de Sitter energy densities 



2 8nGp, 
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Introducing the Planck density ppianck 
the spectrum is given by 



(15) 



hG 2 ' 



1 dp. 



sgu- 



f dp 



sgw 
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PdS 



p c din/ p c df 



9 ppianck 



(16) 



At this point, some comments are in order. First, 
the calculation works for a simplified model that 
does not include the matter-dominated era. If the 
latter is included, the redshift at the equivalence 
epoch has to be considered. Taking into account 
Ref. |4|10J one gets 



^sgw (/) 



16 Pds 

9 PPlanc 



il + Zcg)- 1 



(17) 



for the waves which, at the epoch in which the 
universe becomes matter-dominated, have a fre- 
quency higher than H eq , the Hubble parameter 
at equivalence. This situation corresponds to fre- 
quencies / > (1 + z eq ) 1/2 H a today. The redshift 
correction in eq. (flTI is needed since the present 
value of the Hubble parameter Ho would be dif- 
ferent without a matter-dominated contribution. 
At lower frequencies, the spectrum is given by 



^sgw(f) <X /" 



(18) 



Nevertheless, since the spectrum falls off as f~ 2 
at low frequencies, today at LIGO/VIRGO and 
LISA frequencies, one gets 



^.gw(f)hloo < 2.3 x 1(T 12 



(19) 



where h 100 = Hq/ (100 km • s" 1 • Mpc^ 1 ). It is 
interesting to calculate the corresponding strain 
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at / ~ 100Hz, where interferometers such as 
VIRGO and LI GO achieve maximum sensitivity. 
The well known equation for the characteristic 
amplitude adapted to the scalar component of 
gravitational waves 



$ c (/) ~ 1.26 x 10 



-18 



1H2 

~T 



can be used to obtain 

$ c (100Hz) < 2 - 10" 26 



*J ' h^ 00 H, sgw (/) , 
(20) 

(21) 



Then, since we expect a sensitivity of the or- 
der of 10~ 22 for the above interferometers at 
/ ~ 100 Hz, we need to gain four orders of mag- 
nitude. Let us analyze the situation also at lower 
frequencies. The sensitivity of the VIRGO inter- 
ferometer is of the order of 10 -21 at / ~ 10 Hz 
and in that case it is 



$ c (10Hz) < 2 - 10 



-25 



(22) 



The sensitivity of the LISA interferometer will be 
of the order of 10~ 22 at / ~ 10~ 3 Hz and in this 
case it is 



$ c (l0" 3 Hz) < 2 • 10" 



(23) 



This means that a stochastic background of relic 
scalar gravitational waves could, in principle, be 
detected by the LISA interferometer |llj . 



choose a class of functional forms of f(R) capa- 
ble of matching, in principle, observational data. 
Firstly, the cosmological model should reproduce 
the CMBR constraints in the high-redshift regime 
(which agree with the presence of an effective 
cosmological constant). Secondly, it should give 
rise to an accelerated expansion, at low redshift, 
according to the ACDM model. Thirdly, there 
should be sufficient degrees of freedom in the 
parameterization to encompass low redshift phe- 
nomena (e.g. the large scale structure) accord- 
ing to the observations. Finally, small deviations 
from GR should be consistent with Solar System 
tests. All these requirements suggest that we can 
assume the limits 



lim f(R) = constant, 

R — >oo 



lim f(R) = 0, 



(24) 



(25) 



which are satisfied by a general class of broken 
power law models, proposed in [5], which are 



Jhs(R) — R — Ai? c 



ay 



/ \ 2n 

(A) +1 



(26) 



where m is a mass scale and ci : 2 are dimension- 
less parameters. Besides, another viable class of 
models was proposed in [2] 



5. /(-R)-gravity constrained by PPN pa- 
rameters and stochastic background of 
GWs 

A Brans-Dike-like theory with lo — is di- 
namically equivalent to an /(i?)-gravity, so the 
bounds coming from the interferometric ground- 
based (VIRGO, LIGO) and space (LISA) experi- 
ments could constitute a further probe for gravity 
if matched with bounds at other scales to achieve 
experimental bounds for the theory at local and 
cosmological scales . For our aims consider a class 
of f(R) models which do not contain cosmological 
constant and are explicitly designed to satisfy cos- 
mological and Solar-System constraints in given 
limits of the parameter space. In practice, we 



fs(R)=R + XRc 



1 R 2 

1 + m 



(27) 



Since f(R = 0) = 0, the cosmological constant 
has to disappear in a flat spacetime. The param- 
eters {n, p, A, R c } are constants which should be 
determined by experimental bounds. In Fig.JT]), 
we have plotted the selected models as function 
R 

of — for suitable values of {p, n, A} . 

R c 

The above models can be constrained at Solar 
System level by considering the PPN formalism. 
This approach is extremely important in order to 
test gravitational theories and to compare them 
with GR. One can derive the PPN-parameters 7 



5 




Figure 1. Plots of f(R) models as function of -g- 



Model f HS in Eq. ((26J) with n = 1 and A = 2 
(green line). Model fs in Eq.JSIJ with p = 2, 
A = 0.95 (red line). We also plot f(R)=R (blu 
line) to see whether or not the stability condition 
f r > is violated. 



and /3 in terms of a generic analytic function f(R) 
and its derivative 



7 -l = - 



f"(Rf 



f'(R) + 2/"(i?) 2 ' 



f'(R) ■ f"(R) 



2/'(E) + 3/"(i?) 2 



(28) 



(29) 



dR 



These quantities have to fulfill the constraints 
coming from the Solar System experimental tests. 
They are the perihelion shift of Mercury , the Lu- 
nar Laser Ranging , the upper limits coming from 
the Very Long Baseline Interferometry (VLBI) 
and the results obtained from the Cassini space- 
craft mission in the delay of the radio waves trans- 
mission near the Solar conjunction. 



• Mercury perihelion Shift |27 
1(T 3 

• Lunar Laser Ranging 4/3 
1) x 1CT 3 



0- 1| < 3x 



7 - 3 = (0.7 ± 



• Very Long Baseline Interferometer 
|7- 1| < 4 x 10~ 4 

• Cassini Spacecraft 7-I = (2.1 ±2.3) x 10~ 5 

Take into account the above f(R)- 
models, specifically, we have investigated the val- 
ues or the ranges of parameters in which they 
match the above Solar-System experimental con- 
straints. In other words, we use these models to 
search under what circumstances it is possible to 
significantly address cosmological observations by 
/(i?)-gravity and, simultaneously, evade the local 
tests of gravity [5|. At this point, using the above 
LIGO, VIRGO and LISA upper bounds, cal- 
culated for the characteristic amplitude of GW 
scalar component, let us test the /(-R)-gravity 
models. Taking into account the discussion in 
Sec. we have that the GW scalar component 
is derived considering 



$ = — 



3(7 



CO 

5a = 



a = -ln(l + /'0R))= In f'(A), 



f"(A) 
1 + f'(A) 



5A. 



(30) 



Finally we obtain a good sets of parameters for 
the Starobinsky model PPN vs GW-stochastic: 



• £ = 3.38 ,p=l , A 



• jL = V3 ,p = 2 , 0.944 < A < 0.966 
and for Hu and Sawiki model 



such sets of parameters are not in conflict with 
bounds coming from the cosmological stochas- 
tic background of GWs and, some sets reproduce 
quite well both the PPN upper limits and the 
constraints on the scalar component amplitude 
of GWs. The results indicate that self-consistent 
models could be achieved comparing experimen- 
tal data at very different scales without extrapo- 
lating results obtained only at a given scale [5]. 

6. Conclusions 

We have investigated the possibility that some 
viable f(R) models could be constrained consid- 
ering both Solar System experiments and upper 
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bounds on the stochastic background of gravita- 
tional radiation. Such bounds come from inter- 
ferometric ground-based (VIRGO and LIGO) and 
space (LISA) experiments. The underlying phi- 
losophy is to show that the f(R) approach, in or- 
der to describe consistently the observed universe, 
should be tested at very different scales, that is 
at very different redshifts. In other words, such a 
proposal could partially contribute to remove the 
unpleasant degeneracy affecting the wide class of 
dark energy models, today on the ground. 

Beside the request to evade the Solar System 
tests, new methods have been recently proposed 
to investigate the evolution and the power spec- 
trum of cosmological perturbations in f(R) mod- 
els. The investigation of stochastic background, 
in particular of the scalar component of GWs 
coming from the f(R) additional degrees of free- 
dom, could acquire, if revealed by the running 
and forthcoming experiments, a fundamental im- 
portance to discriminate among the various grav- 
ity theories [9|. These data (today only up- 
per bounds coming from simulations) if combined 
with Solar System tests, CMBR anisotropics, 
LSS, etc. could greatly help to achieve a self- 
consistent cosmology bypassing the shortcomings 
of ACDM model. 

Specifically, we have taken into account two 
broken power law f(R) models fulfilling the main 
cosmological requirements which are to match 
the today observed accelerated expansion and the 
correct behavior in early epochs. We have taken 
into account the results of the main Solar Sys- 
tem current experiments. Such results give up- 
per limits on the PPN parameters which any self- 
consistent theory of gravity should satisfy at local 
scales. Starting from these, we have selected the 
f(R) parameters fulfilling the tests. As a general 
remark, all the functional forms chosen for f(R) 
present sets of parameters capable of matching 
the two main PPN quantities, that is j exp and 
flexp- This means that, in principle, extensions of 
GR are not a priori excluded as reasonable can- 
didates for gravity theories. To construct such 
extensions, the reconstruction method developed 
in may be applied. 

The interesting feature, and the main result of 
this paper, is that such sets of parameters are not 



in conflict with bounds coming from the cosmo- 
logical stochastic background of GWs. In particu- 
lar, some sets of parameters reproduce quite well 
both the PPN upper limits and the constraints on 
the scalar component amplitude of GWs |5|llj . 

Far to be definitive, these preliminary re- 
sults indicate that self-consistent models could 
be achieved comparing experimental data at very 
different scales without extrapolating results ob- 
tained only at a given scale. 
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Table 1 



The next-to-leading 


; order (NLO) results 


without the 


pion field. 




A (McV) 


140 


150 


175 


200 


r d (fm) 
Qd (fm 2 ) 
Pd (%) 
f'd 

Mmi (fm) 
A^gt (fm) 

^ (%) ' 

S™> (%) 
S™" (%) 


1.973 

0.259 

2.32 

0.867 

3.995 

4.887 

-0.45 
0.03 

-0.19 


1.972 

0.268 

2.83 

0.864 

3.989 

4.881 

-0.45 
0.03 

-0.19 


1.974 

0.287 

4.34 

0.855 

3.973 

4.864 

-0.45 
0.03 

-0.18 


1.978 

0.302 

6.14 

0.845 

3.955 

4.846 

-0.45 
0.03 

-0.15 



The experimental values are given in ref. [2]. 



packages can be found in I^TgX's tools directory. 

array Various extensions to I^Tj^X's array and 
tabular environments. 

longtable Automatically break tables over sev- 
eral pages. Put the table in the longtable 
environment instead of the table environ- 
ment. 

dcolumn Define your own type of column. 
Among others, this is one way to obtain 
alignment on the decimal point. 

tabularx Smart column width calculation 
within a specified table width. 

rotating Print a page with a wide table or 
figure in landscape orientation using the 
sidewaystable or sidewaysf igure envi- 
ronments, and many other rotating tricks. 
Use the package with the f iguresright op- 
tion to make all tables and figures rotate 
in clockwise. Use the starred form of the 
sideways environments to obtain full-width 
tables or figures in a two-column article. 

3.3. Line drawings 

Line drawings may consist of laser-printed 
graphics or professionally drawn figures attached 
to the manuscript page. All figures should be 
clearly displayed by leaving at least one line of 
spacing above and below them. When placing a 



figure at the top of a page, the top of the figure 
should align with the bottom of the first text line 
of the other column. 

Do not use too light or too dark shading in your 
figures; too dark a shading may become too dense 
while a very light shading made of tiny points may 
fade away during reproduction. 

All notations and lettering should be no less 
than 2 mm high. The use of heavy black, bold 
lettering should be avoided as this will look un- 
pleasantly dark when printed. 



3.4. PostScript figures 

Instead of providing separate drawings or 
prints of the figures you may also use Post- 
Script files which are included into your DTj^X 
file and printed together with the text. Use 
one of the packages from DT^X's graphics di- 
rectory: graphics, graphicx or epsfig, with 
the \usepackage command, and then use the 
appropriate commands (\includegraphics or 
\epsf ig) to include your PostScript file. 

The simplest command is: 
\includegraphics{f ile}, which inserts the 
PostScript file file at its own size. The starred 
version of this command: 

\includegraphics*{f ile}, does the same, but 
clips the figure to its bounding box. 

With the graphicx package one may spec- 
ify a series of options as a key-value list, e.g.: 
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\includegraphics [width=15pc] {file} 
\includegraphics [height=5pc] {file} 
\includegraphics [scale=0 . 6] {file} 
\includegraphics [angle=90,width=20pc] {file} 

See the file grfguide, section "Including 
Graphics Files", of the graphics distribution for 
all options and a detailed description. 

The epsf ig package mimicks the commands fa- 
miliar from the package with the same name in 
MeX2.09. A PostScript file file is included with 
the command \psf ig{f ile=f ile}. 

Grey-scale and colour photographs cannot be 
included in this way, since reproduction from the 
printed CRC article would give insufficient typo- 
graphical quality. See the following subsections. 



Figure 2. Remember to keep details clear and 
large enough. 



Figure 1. Good sharp prints should be used and 
not (distorted) photocopies. 



3.5. Black and white photographs 

Photographs must always be sharp originals 
(not screened versions) and rich in contrast. 
They will undergo the same reduction as the text 
and should be pasted on your page in the same 
way as line drawings. 

3.6. Colour photographs 

Sharp originals (not transparencies or slides) 
should be submitted close to the size expected 
in publication. Charges for the processing and 
printing of colour will be passed on to the 



author(s) of the paper. As costs involved are per 
page, care should be taken in the selection of size 
and shape so that two or more illustrations may 
be fitted together on one page. Please contact 
the Author Support Department at Elsevier (E- 
mail: authorsupport@elsevier.nl) for a price 
quotation and layout instructions before produc- 
ing your paper in its final form. 

4. EQUATIONS 

Equations should be flush-left with the text 
margin; F/IgX ensures that the equation is pre- 
ceded and followed by one line of white space. 
ET^X provides the document class option f leqn 
to get the flush-left effect. 

= EW(u)5 af3 + (a\W n \(3) (1) 

You need not put in equation numbers, since 
this is taken care of automatically. The equation 
numbers are always consecutive and are printed 
in parentheses flush with the right-hand margin 
of the text and level with the last line of the equa- 
tion. For multi-line equations, use the eqnarray 
environment. 

For complex mathematics, use the A\4Sma,th 
package. This package sets the math indenta- 
tion to a positive value. To keep the equations 
flush left, either load the espcrc package af- 
ter the .A^Smath package or set the command 



\mathindent=Opt in the preamble of your arti- 
cle. 
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References should be collected at the end of 
your paper. Do not begin them on a new page 
unless this is absolutely necessary. They should 
be prepared according to the sequential numeric 
system making sure that all material mentioned 
is generally available to the reader. Use \cite 
to refer to the entries in the bibliography so that 
your accumulated list corresponds to the citations 
made in the text body. 

Above we have listed some references according 
to the sequential numeric system |1I2I3I4| . 



